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Rocket Propulsion: Volcanic Caldera on Titan 
(A Simulation) 


1 Introduction 


In this application note, the mission is to determine the speed of a rocKet- 
powered drone sled, as it accelerates across a volcanic caldera icy surface. 
Quadratic drag is assumed to be the only opposing force while sliding friction 
is ignored. 


Initially, the sled is at rest. After the engine is ignited (powered phase), the 
sled’s speed increases until engine cut off at 20 seconds. Then, the sled simply 
coasts. However, the coasting phase will not be addressed in this simulation. 


For comparison, analytical, numerical, and analog computation solutions are 
displayed in table 1. 


2 Mathematical modeling 


Rocket propulsion equation derivation: 
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Figure 1: Rocket Propulsion 


where, 
m= mass of rocket and fuel 
—Am = incremental mass of ejected fuel (note that Am < 0) 
v =speed of sled relative to a fixed frame of reference 
Av = incremental speed increase of the rocket and remaining fuel 
v = speed of ejected fuel relative to a fixed frame of reference 
u = speed of ejected fuel relative to the sled (assumed constant) 
Fext = external forces on rocket (in this case, quadratic drag) 
Starting with impulse of the force = change in momentum, 
FextAt = Ap 
FextAt = Patter — Pbefore 
FextAt = [(m + Am)(v + Av) + (-Am)(v’)] - [mv] 
FextAt = [(m + Am)(v + Av) + (-Am)(v - u)] - [mv] 


FexAt = [mv + mAv + vAm + AmAv - vAm + uAm | - [mv] 





FexrAt = mAv + AmAv + uAm 





Fext = mAv/At + AmAv/At + uAm/At 





Fext = mAv/At + (Am/At)(Av/At) At + uAm/At 


lim (Fext) = Jim [mAv/At + (Am/At)( Av/At) At + u Am/At] 
lim (Fext) = mdv/dt + udm/dt 


Assuming Fext = -quadratic drag = -0.5CppAv?, 


jim (—0.5CppAv*) = mdv/dt + udm/dt 


—~0.5CppAv2 = mdv/dt + udm/dt 


After a little algebra, 
dv/dt = (—0.5CppAv? — udm/dt)/m 
Letting dm/dt = —c (constant), where = c > 0 and m(0) = mo, 
m= -—ct + Mo 
Assuming u is constant, 
dv/dt = (—0.5CppAv?2 + uc)/(—ct + mo) 


If mo = 1.800 kg, c = 0.075, kg/s, Co= 0.250, p = 4.000 kg/m3, A = 0.100 m2, 
u = 24 m/s, and omitting units for simplicity, then 


dv/dt = (-0.050v?2 + 1.800)/(—0.075t + 1.800) with v(0) = 0 m/s 


Use the following equation for the analog computer setup and numerical 
methods: 


dv/dt = (—0.050v? + 1.800) /(-0.075t + 1.800) 


For mathematical convenience to avoid decimals, multiply the numerator and 
denominator of the above equation by -20. 


Simplifying, 
dv/dt = (v2 — 36)/(3/2 t — 36) 
Multiplying the numerator and denominator by 2/3, 
dv/dt = 2/3(v2- 62)/(t— 24) (1) 
Setting dv/dt = 0 to determine the sled’s terminal speed, 


0 = 2/3(v4— 62)/(t— 24) 


By inspection, vr = 6. 
Returning to (1), 
dv/dt = 2/3(v? — 62)/(t — 12), which is separable and integrable. 


dv /(v2 — 62) = 2/3 dt/(t— 24) 


dv/(v? — 62) = 2/3 i (2) 
/< 


Using the method of partial fractions for the integrand on the left, 
1/(v2 - 62) = A/(v- 6) + B/(v + 6) (3) 
OV+1=(A+B)v + 6(A -B) 

Comparing coefficients of like terms, 
A+B=0O0orB=-A 
and 
6(A - B) = 6(A+ A) = 1 0rA=1/12 and B=-1/12 
Substituting the values for A and B into (3), 

1/(v2 - 62) = 1/12[1/(v - 6) - 1/(v + 6)]| (4) 


Substituting (3) into the integrand on the left side of (1), 


Vv t 


1/12 Jae ~6) —1/(v + 6)]dv = 2/3 i dt/(t — 24) (5) 


0 


Jive ~6)—-1/(v + 6)]dv =8 | dt/(t — 24) (6) 


Integrating (6) by inspection, 
In|(v- 6)/(v + 6)||9 = 8 In|t - 24|]6 
In|(v - 6)/(v+6)| = 8 In]1 - t/24| 
|(v - 6)/(v+6)| = e*(8 In|1 - t/24]) 
Since v < 6 andt < 20, 
—(v - 6)/(v+6) = e*(8 In(1 - t/24)) 
Letting E = e*(8 In(1 - t/24)) for simplicity, 
—(v - 6)/(v+6) =E 
After some algebra, 
v=6(1-E)/(1 + E) 
Since E > 0, then (1 - E)/(1 + E) < 1. Confirming that v < 6. 
Replacing e*(8 In(1 - t/24) for E, 


3a Computer setup (patch cord version) 
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Figure 2: Computer setup for Rocket Propulsion: Volcanic Caldera on Titan 





3b Computer setup (op amp/discrete component version) 


to linear function generator 


voltage divider 


to linear function generator 
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Arrow Color Code: 


Red arrows connect to Vcc. 
Blue arrows connect to Vee. 
Green arrows connect to ground. 


Resistance values in kQ 
Subcircuit Color Code: 
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Inverting Integrator 


Red 
Analog multiplier 
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Inverting Summer 
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Figure 3: Basic Breadboard Layout 


4 Numerical Method 
(Modified Euler method using a hand-held programmable calculator) 
Code: TI-BASIC 


PROGRAM:CALDERA 
:-ClrHome:ClrDraw 
’DV/DT=(-0.050V2+UC) /(-CT+M)” Note: M = initial mass 
“INITIAL CONDITION:” 
:(0-T:0-V 

"ENGINE PARAMETERS:” 
:24-U:0.075-C 
:1.800-M:20-E 

:”’STEP SIZE:” 

:0.100-—H 

Fix 2 

:Lbl 1 

‘If T>E:Then 

:Goto 2:Else 

:Disp {T,V} 
:(-0.050V2+UC)/(-CT+M)>F 
‘T+HOT 

:V+FH>W 

: (-0.050W2+UC)/(-CT+M)-S 
((F+S)/2-A 

:V+AH>V 

:Pause 

:Gotol 

:Lbl 2 

‘End 


5 Results 



































t(s) Analog computer Analytical Numerical Method 
v(m/s) v(m/s) v(m/s) 
estimated from Calculator 
oscilloscope 

00.0 0.00 0.00 0.00 

02.5 2.48 2.48 2.48 

05.0 4.40 4.40 4.40 

07.5 5.40 5.43 5.43 

10.0 5.90 5.84 5.84 

125 6.00 5.97 5.97 

15.0 6.00 6.00 6.00 

17.5 6.00 6.00 6.00 

20.0 6.00 6.00 6.00 

















Table 1: Solution comparisons 





Figure 4: Drone sled speed on Titan* 


*For this application note, the oscilloscope display was produced during a 
single run using a differential equation analog computer constructed from 
operational amplifiers and discrete components with tolerances within 10%. 
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Figure 4: Differential Equation Analog Computer** 


** The photograph does not include the DC supply with connections. 
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